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Abstract 

We analyze the four-fermion reactions e + e~ — * 4/ containing a single top quark and three 
other fermions, a possible decay product of the resonant anti-top quark, in the final state. 
This allows us to estimate the contribution of the nonresonant Feynman graphs and effects 
related to the off mass shell production and decay of the top quark. We test the sensitivity 
of the total cross section at centre of mass energies in the ti threshold region and far above it 
to the variation of the top quark width. We perform calculation in an arbitrary linear gauge 
in the framework of the Standard Model and discuss an important issue of gauge symmetry 
violation by the constant top quark width. 
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1 Introduction 



Production of the top quark in e + e~ annihilation is an issue which has attracted a lot of interest 
of both experimenters and theorists for the past decade. In the present decade the interest will 
be certainly growing up in the prospect of new high energy e + e~ colliders gradually coming 
up from the stage of general discussion to a detailed planning and hopefully to the stage 
of construction and successful operation [|l|. As the e + e~ machines operate in a very clean 
experimental environment they provide a unique possibility of precise measurement of the top 
quark physical properties which, as many expect, may go beyond the standard model (SM) 
and give decisive hints for the development of new physical ideas. In order to disentangle 
the possible effects of new physics from the standard physics it is crucial to know the SM 
predictions as precisely as possible. 

It has taken a lot of effort to obtain precise SM predictions for the top pair production in 
the threshold region. A substantial improvement of convergence of the perturbation series has 
been achieved by computing the next-to-next-to-leading order corrections to the top quark 
pair production cross section near threshold || and understanding the renormalon cancellation 
mechanism H. The 0(aa s ) corrections to the top decay with into a W boson and a b quark 
are also known ||]. 

In the present note, we concentrate on the effects related to the fact that one of the quarks 
in the ti pair may be produced off mass shell. In particular, we let the anti-top quark i decay 
into a final state possible in the framework of the SM, i.e., we consider reactions 

e+e" - tiff, (1) 

where / = e~ , \T~ ,r~ ,d, s and /' = u e , z/ M , u T , u, c, respectively, taking into account the com- 
plete set of the Feynman graphs which contribute to the specific final state at the tree level. 
We pay attention to the very important issue of gauge symmetry breaking caused by the 
nonzero constant widths of unstable particles, in particular that of the top quark, which are 
kept as free parameters. We also test the sensitivity of the total cross sections of ([I]) to the 
variation of the top width. The deviation of the top width from its SM value may indicate the 
new physics. In order to test the reliability of our results, we perform the calculation in the 
arbitrary linear gauge. We neglect radiative corrections, the correct treatment of which 
demands an extra effort and is beyond the scope of the present work. 

We describe basics of the calculation in the next section. Our results are presented and 
discussed in Section 3 and, finally, in Section 4, we give our concluding remarks. 



2 Calculation in arbitrary linear gauge 

The calculation of the necessary matrix elements relies on the method proposed in Ref. || 
and further developed in Ref. 0. As in Ref. ||, fermion masses are kept nonzero both in the 
matrix elements and in the kinematics. Among others, this has an advantage that the Higgs 
boson effects can be incorporated consistently and a pole related to the photon exchange in 
the t-channel can be handled better than in the massless fermion case. 

In order to estimate the gauge symmetry violation effects related to the nonzero widths of 



unstable particles, we perform the calculation in two different schemes: the 'fixed widths 
scheme' (FWS) and in the so called 'complex-mass scheme' (CMS) of Ref. |7]]. Both schemes 
introduce the constant particle widths through the complex mass parameters: 



M^ = m 2 v - im v T v , V = W,Z, H, and M t = m t - iT t /2, (2) 

which replace masses in the corresponding propagators. The coupling constants are given in 
terms of the electric charge and electroweak mixing parameter sin 2 #iy. In FWS, the elec- 
troweak mixing parameter is kept real, i.e., it is given by 

sm 2 6 w = 1 - rn%r/m%, (3) 

where mw and mz are physical masses of the W ± and Z° boson. In CMS, sin 2 Ow is given in 
terms of the complex masses Mw and Mz of Eq. (fj) by 

sm 2 6 w = l-M^/M 2 z , (4) 

thus it is a complex number. The CMS has the advantage that it preserves the Ward identities 
0, provided all the fermion widths are zero. We have tested the gauge invariance numerically 
for all the reactions (|l|) under the assumption of the zero top quark widths. With the nonzero 
top width introduced in the Feynman propagator of the top quark 



i8f{p) = i^h, (5) 



p 2 - Mi 

the Ward identities are not satisfied any more and hence the gauge symmetry is violated. 



There has been a lot of discussion of the gauge invariance issue in literature of the past few years 
||. The best way to solve the problem, however involved and complicated it might be, is to 
work in the 'fermion-loop scheme' (FLS), which resumes higher order effects coming from one- 
loop fermion contributions to bosonic propagators together with parts of the vertices necessary 
for keeping the corrections gauge invariant. Unfortunately, a similar Dyson resummation has 
not been worked out in a detailed way for a propagator of an unstable fermion up to now. 
Therefore, the substitution of Eq. @ has no full theoretical justification in the framework of 
quantum field theory. However, a comparison of numerical results obtained in the schemes 
using fixed widths in both the s- and t-channel gauge boson propagators with those derived 
within FLS, which show no numerically relevant discrepancies for several four-fermion reactions 
and the corresponding bremsstrahlung processes, speaks in favour of the simplified approach 
based on substitution of Eq. (^|) in each propagator of unstable particle. Therefore, in the 
following we will restrict ourselves to this simplified approach. 

In order to quantitatively estimate gauge symmetry violation effects induced by the substi- 
tution of Eq. (FJ), we perform calculation of the matrix elements in the linear gauge with 
arbitrary real gauge parameters £y, V = j,W,Z. We then vary the parameters in a very wide 
range from £y = 1 corresponding to the 'tHooft-Feynman gauge (FG) to £y = 10 16 which, 
in the double precision of Fortran programming language, corresponds to the unitary gauge 
(UG). We would like to stress at this point that we take into account contributions from the 
exchange of the would-be Goldstone bosons in the gauge, which are absent in the unitary 
gauge. If the change in the cross section induced by the change in gauge parameters is smaller 
than the accuracy of the Monte Carlo integration we assume that the gauge violation effects 
induced by the nonzero top width are numerically irrelevant and we may consider the corre- 
sponding results as trustworthy. On the other hand, if the results depend on the choice of 



gauge parameters they are useless, but we may try to reduce the dependence by imposing cuts 
on the phase space integration. This simple prescription, however doubtful from the purely 
theoretical point of view might it be, allows one to treat the particle widths as independent 
parameters and test the reliability of the numerical results in an efficient way. 

The phase space integration is performed numerically using a multichannel Monte Carlo (M C) 
approach and the integration routine VEGAS [|l(J. The 7 dimensional phase space element of 
reaction ([]]) is parametrized in a few different ways in order to account for the most relevant 
peaks of the matrix elements: the ~ 1/t pole caused by the t-channel photon-exchange, the 
Breit-Wigner shape of the W and Z° resonances, the ~ 1/s behavior of a light fermion pair 
production, and the ~ 1/t pole due to the the neutrino exchange at the same time. 



3 Numerical results 

In this section, we will present numerical results for all the four-fermion channels of reaction 
([l|) possible in the SM. 

We define the SM physical parameters in terms of the gauge boson masses and widths, the 
top mass and the Fermi coupling constant. We take the actual values of the parameters from 
Ref. §: ' 

m w = 80.419 GeV, T w = 2.12 GeV, m z = 91.1882 GeV, Y z = 2.4952 GeV, 

m t = 174.3 GeV, = 1.16639 x 1(T 5 GeV" 2 . (6) 

We assume the Higgs boson mass of mu = 115 GeV and, if not stated otherwise, the top 
quark width of T t = 1.5 GeV. 

For the sake of definiteness we also list other fermion masses we use in the calculation ||: 

m e = 0.510998902 MeV, = 105.658357 MeV, m T = 1777.03 MeV, 

m u = 5 MeV, m d = 9 MeV, m s = 150 MeV, m c = 1.3 GeV, m b = 4.4 GeV. (7) 

We neglect the Cabibo-Kobayashi-Maskawa mixing, i.e., we assume the CKM matrix to be a 
unit matrix. 

The fine structure constant is calculated from 

a w = \p2G Pu m\ v sin 2 6 W / n (8) 
with the real electroweak mixing parameters of Eq. (|||) in the both schemes FWS and CMS. 

Except for the check of gauge invariance discussed in the previous section, we perform a few 
other checks. Our results reproduce those of Ref. for a top mass smaller than my/ and 
the zero top width. The corresponding matrix elements in the absence of the Higgs boson 
exchange has been checked against MADGRAPH [[□]]. The phase space generation routine for 
particles of large masses has been written in two independent ways. 

In Table 1, we show the results for the cross sections of e + e~ — > tb^v^ at different centre 
of mass energies obtained in different schemes and gauges: the complex-mass scheme (CMS), 



fixed width scheme (FWS), unitary gauge (UG) and Feynman gauge (FG). We have integrated 
over the full four particle phase space without any cuts. We can see that the results hardly 
depend on the gauge choice both in the CMS and FWS. Actually, they nicely agree with each 
other within one standard deviation of the MC integration. 

In Table 2, we present the results for e + e~ — > tbe~u e obtained in the CMS and in two different 
gauges, UG and FG. In order to reduce the dependence on gauge choice induced by the nonzero 
top width, we have imposed a cut on the electron angle with respect to the beam 6(e~, beam). 
Again there is rather small dependence on the cut for the energies presented in Table 2. From 
a comparison with the corresponding numbers of Table 1, we can infer that the t-channel 
Feynman graphs of reaction e + e~ — > tbe~u e do not contribute much to the total cross section 
in the presence of the cut on the final electron angle. When we reduce the cut further so that 
the denominator of the t-channel photon propagator becomes of the order of the electron mass 
squared the dependence on gauge becomes substantial and results are meaningless. 



Table 1: Cross sections in fb of e + e~ — > tbyTv^ at different centre of mass energies in 
different schemes, CMS and FWS, and gauges, UG and FG. The numbers in parenthesis 

show the uncertainty of the last decimals. 



(GeV) 


UG 
a CMS 


FG 
a CMS 


UG 
a FWS 


CT FG 

a FWS 


190 


2.6174(7) x 10~ 8 


2.6174(7) x 10~ 8 


2.6185(7) x 10~ 8 


2.6185(7) x 10~ 8 


340 


0.7837(4) 


0.7837(4) 


0.7839(4) 


0.7840(4) 


360 


41.27(10) 


41.27(10) 


41.28(10) 


41.29(10) 


500 


60.06(13) 


60.04(13) 


59.75(30) 


59.90(29) 


2000 


5.59(3) 


5.56(3) 


5.51(7) 


5.51(8) 



Table 2: Cross sections in fb of e + e — > the u e in the CMS in two different gauges, UG and 
FG, and for two different cuts on the electron angle with respect to the beam. 



(GeV) 


5° < #(e-,beam) < 175° 


1° < fl(e-,beam) < 179° 


UG 
a CMS 


CT FG 

a CMS 


UG 
a CMS 


CT FG 

a CMS 


190 
340 
360 
500 
2000 


0.6607(4) x 10~ 5 
0.7993(4) 
41.21(11) 
59.78(15) 
6.81(3) 


0.6607(4) x 10~ 5 
0.7993(4) 
41.20(11) 
59.75(15) 
6.82(3) 


0.10520(4) x 10~ 4 
0.8251(4) 
41.32(8) 
60.16(15) 
7.97(3) 


0.10531(4) x 10" 4 
0.8253(4) 
41.32(8) 
60.19(15) 
8.00(3) 



The results for the channels of reaction ([I]) which do not contain electron in the final state are 
shown in Table 3. They were obtained in the CMS and unitary gauge. 



Table 3: Cross sections in the CMS in fb of different channels of reaction (|T|) not containing 

a final state electron. 



Channel of 
reaction (|l|) 


yfi (GeV) 


190 


340 


360 


500 


2000 


e + e~ — > tbfi'u^ 


2.6174(7) x 1(T 8 


0.7837(4) 


41.3(1) 


59.8(3) 


5.42(7) 


e + e~ — > tbr"v T 


1.9331(4) x 1(T 8 


0.7831(4) 


41.2(1) 


59.6(3) 


5.47(7) 


e + e~ — > tbdu 


7.880(2) x 10- 8 


2.351(1) 


123.8(3) 


179.9(9) 


16.3(2) 


e + e~ — > tbsc 


6.616(2) x 10~ 8 


2.350(1) 


123.8(3) 


178.9(9) 


16.7(2) 



In Fig. 1, we show the energy dependence of the total cross section of e + e~ — ► tbyTv^ calculated 
with the complete set of Feynman graphs and the approximate cross section e + e~ — > ti — > 
tbyTVu. The latter has been obtained by multiplying the on shell top pair production cross 
section by the corresponding three body top decay width 

a(e + e~ — > ti — > tbfi^Ufj) = cr(e + e~ — > ti) T(i — > b]i~u^). (9) 

We have taken over the SM part of the analytic formula for the width T(i —>■ bpTv^) with 
massless final state fermions from Ref. |l!|. In the calculation of a(e + e~ — > tb/jTv^) we have 
used the FWS scheme and neglected the Higgs boson contribution. We see that Eq. ([]) 
approximates the complete tree level calculation well not only just above the threshold, but 
also for higher centre of mass energies. The relative difference between the both results is 
3.5% at 360 GeV, 1.3% at 500 GeV and - 5.0% at 800 GeV. The nice agreement is somewhat 
amazing as except for one Feynman graph which contain a resonant top propagator there are 
nine other nonresonant graphs which contribute to e + e~ — > thyTv^ in the unitary gauge. 

The explanation of this fact can easily be found if one looks at Fig. 2 where we have plotted, 
against the centre of mass energy, the cross section of Eq. (|9]) and another approximated cross 
section obtained by integrating over full four particle phase space the squared matrix element 
containing only the top resonant Feynman graph. The small discrepancy between the two 
curves in Fig. 2 is a measure of spin correlations and off shellness of the i quark. 

We illustrate the dependence of the total cross section of e + e~ — > thyTv^ on the top quark 
width r t in Table 4. We see that the cross section at ^fs = 360 GeV, i.e. just above the 
threshold is almost exactly proportional to 1/IY This kind of dependence holds also at 
y/s = 500 GeV, which is already much above the threshold and it survives almost unaltered 
at y/s = 2 TeV. It means that the cross section of e + e~ — > ibjjrv^ is well approximated by the 
resonant i production and its subsequent decay. This kind of dependence offers a new way of 
measurement of the top quark width alternative to the measurement based on the shape of 
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Figure 1. The energy dependence of the total cross sections of e + e 

e + e~ — ► ti — > tbfi~Un. 
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Figure 2. The energy dependence of the total cross sections of e + e~ - 

e + e~ — ► ti — ^ tb^L~v^. 



Table 4: Cross sections in fb of e + e — ► tb/i for different values of the top quark width. 
The calculation has been performed in the CMS and UG. 



(GeV) 


r t (GeV) 


1.5 


1.6 


1.7 


190 


2.6174(7) x 10- 8 


2.6186(4) x 10" 8 


2.6186(4) x 10" 8 


340 


0.7837(4) 


0.7832(3) 


0.7830(3) 


360 


41.27(10) 


38.65(7) 


36.31(6) 


500 


60.06(13) 


56.37(13) 


53.13(12) 


2000 


5.59(3) 


5.30(2) 


5.09(2) 



4 Summary and Outlook 

We have analyzed the top quark production in e + e~ annihilation at a new high luminosity 
linear collider like TESLA. We have estimated the contribution of the nonresonant Feynman 
graphs and effects related to the off mass shell production and decay of one of the top quarks. 
Those effects are typically of the order of a few per cent. Therefore one should take them into 
account in the analysis of the future data. We have shown that the cross section of reaction 
([]]) is dominated by the resonant i production and its subsequent decay not only at centre of 
mass energies in the tt threshold region but also far above it. We have tested the sensitivity 
of the total cross sections to the variation of the top quark width and confirmed expected 
proportionality to 1/T t over a very wide energy range beginning from the threshold. This 
kind of dependence offers an alternative way of measurement of the top quark width. By 
performing calculation in an arbitrary linear gauge in the framework of the Standard Model 
we have been able to address an important issue of gauge symmetry violation by the constant 
top quark width. 

It would be desirable to consider the effects related to the off shell production of the second 
quark of the tt pair and to include leading radiative corrections in the analysis of the future 
data. 
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